UK  UNLIMITED 


/ 


Procurement  Executive  -  Ministry  of  Defence 


Recommended  for  issue  by 

P  P  A  Richards,  Superintendent 
Approved  by 

F  E  Whiteway,  Head  of  Division 


OQ-U  O  I 

UK  UNLIMITED 


(c)  Controller,  Her  Majesty '.s  Stationery  Office,  London  -  1  980 


ISBN  0  855181^1  9 


CONTENTS 


SUMMARY 

INTRODUCTION 

THE  BASIC  GENERATOR  MODEL 

Generator  inductance 
Armature  effective  resistance 

Equivalent  skin  depth  method 
Asymptotic  analytic  method 
Finite  difference  method 

Coil  effective  resistance 

BASIC  INPUTS  TO  THE  MODEL 

Conductor  properties 
Armature  dynamics 

CODE  RESULTS 

Typical  results 
Compatibility  of  models 
Load  effects 

Influence  of  wire  diameter  and  winding  pitch 
Influence  of  expansion  angle  and  expansion  ratio 
Influence  of  priming  current 
Influence  of  resistivity 

CONCLUSIONS 

REFERENCES 

APPENDIX  A 

FIGURES  1  -  22 


Accession  For 

lirii  gHU&i 

MX)  TAB 

Unatmounced 

Justification 


Distrihuti-^r/ _ 

Aval  lab  j ”  d 

Aval  1  -d/ci' 
Dlat.  special 


SUMMARY 


An  equivalent  circuit  model  is  presented  for  a  helical  compressed 
magnetic  field  generator.  The  emphasis  has  been  placed  on  producing  a  model 
which  has  a  short  computer  run  time. 

Magnetic  energy  losses  due  to  non-linear  field  diffusion  are  taken  into 
account,  together  with  the  effects  of  magnetically-induced  conductor  motion. 

Examples  of  the  computer  code  results  are  given,  together  with  a 
comparison  with  experimental  data. 


1.  INTRODUCTION 

The  amplification  of  magnetic  fields  due  to  the  driving  of  current 
carrying  conductors  at  high  velocities  has  found  a  number  of  applications.  For 
instance,  ‘he  high  energy  densities  obtainable  have  been  used  to  heat  plasmas. 
Alternatively,  the  corresponding  amplification  of  the  current  in  the  associated 
conductors  serves  as  the  basis  for  electrical  energy  generation  in  low  impedance 
loads.  The  author  is  mainly  concerned  with  this  latter  application. 

In  order  to  transduce  significant  quantities  of  electrical  energy  large 
conductor  velocities  are  necessary.  Too  small  a  velocity  leads  to  large  diffusion 
losses  resulting  in  a  poor  or  even  to  an  absence  of  energy  gain.  Conductors  can  be 
accelerated  to  substantial  velocities  by  subjecting  them  to  explosive  impulse. 
Modern  explosive  technology  permits  controlled  and  repeatable  experimentation 
in  this  area. 

The  helical  generator  allows  one  to  conveniently  match  the  device  to 
the  load  and  is  the  class  of  generator  under  consideration  here.  Such  a  device  is 
shown  schematically  in  figure  1.  It  consists  of  a  cylindrical  conducting  shell 
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(armature)  containing  a  high  explosive  together  with  a  concentric  helical  coil. 
The  generator  is  primed  (given  an  initial  current/field)  by  discharging  a  capacitor 
into  the  coil.  Synchronised  initiation  of  the  explosive  charge  at  one  end  leads  to  a 
funnelling  of  the  armature  and  the  trapping  of  magnetic  flux  between  the 
armature  and  coil.  The  expanding  armature  systematically  wipes  out  the  space 
between  the  armature  and  coil  leading  to  field,  current  and  magnetic  energy 
amplification.  With  suitable  generator  design,  the  transduced  energy  can  be 
efficiently  coupled  into  the  load. 

Because  of  its  dependence  on  explosive  technology  the  generator  is  of 
a  "one  shot"  nature.  However,  this  disadvantage  is  offset  by  the  generator's 
ability  to  transduce  substantial  quantities  of  electrical  energy  at  small  volume. 
Nevertheless,  because  of  the  destructive  nature  of  the  experiments,  generator 
development  can  be  time-consuming  and  expensive.  A  realistic  model  represents 
an  invaluable  tool  for  minimising  the  time  and  cost  of  development.  The  model 
should  ideally  be  capable  of  identifying  the  significance  of  the  generator's 
parameters  with  a  view  to  carrying  out  optimisation  and  sensitivity  analyses.  As 
a  result,  the  experimental  programme  is  minimised. 

In  order  to  be  effective,  the  model  should  be  capable  of  dealing  with 
the  essential  physics  of  the  generator.  However,  too  rigorous  a  mathematical 
treatment  of  this  coupled  magneto-hydrodynamic  problem  leads  to  very  lengthy 
computer  calculations.  Such  treatments  have  been  known  to  result  in  calculation 
times  of  many  hours.  Although  giving  greater  exactness,  these  models  tend  to  be 
of  more  limited  use  in  parametric  studies  which  entail  large  numbers  of 
calculations.  The  approach  presented  here  aims  at  producing  a  model  containing 
all  the  essential  physics  while  still  maintaining  a  reasonably  fast  computer 
capability. 


The  two  output  characteristics  of  the  generator  which  are  of  most 

interest  are  its  amplification  factors;- 

peak  generated  current 
Current  gain  =  initial  current 

Energy  gain  =  Peak  generated  load  magnetic  energy 

initial  magnetic  energy 


The  physical  generator  design  may  vary  from  that  shown  sche¬ 
matically  in  figure  i.  For  instance,  in  order  to  shorten  the  pulse  shape  the 
explosive  charge  may  be  initiated  simultaneously  at  both  ends  (armature  funnels 
in  from  both  ends)  or  initiated  at  the  mid-point  (armature  funnels  towards  both 
ends).  The  shorter'.ing  of  the  running  time  does  not  lead  to  an  enhancement  of  the 
amplification  factors. 

There  are  two  distinct  phases  in  the  generator's  operation.  The  initial 
phase  involves  the  conversion  of  explosive  energy  into  kinetic  energy  of  the 
armature.  This  kinetic  energy  is  then  transduced  into  magnetic  energy  by  virtue 
of  the  current  carrying  conductor  doing  work  on  the  electromagnetic  forces 
present.  The  high  energy  density  of  explosives  represents  the  key  to  the 
generator's  potential.  However,  the  efficiency  of  conversion  of  kinetic  to 
magnetic  energy  decreases  with  decreasing  size  of  the  generator  because  of  the 
more  severe  diffusion  and  joule  heating  losses  associated  with  small  generators. 
In  larger  systems  high  efficiencies  resulting  from  a  total  transfer  of  armature 
kinetic  energy  (turnaround)  are  achievable. 

The  model  described  in  this  report  deals  with  fixed  diameter  (wire) 
round  wire  constant  pitch  helical  coils  only.  In  general  as  generators  are  scaled 
up  in  size  and  are  run  at  higher  energy  levels,  coils  are  graded  in  order  to  reduce 
losses  due  to  joule  heating  and  magnetic  pressure  effects  (ie,  wire  diameter  and 
pitch  increase  towards  output  end). 

A  previous  desc-iption  of  this  model  has  been  given  in  reference  (l). 


2.  THE  BASIC  GENERATOR  MODEL 

This  takes  the  form  of  a  lumped  equivalent  electrical  circuit  and  is 
shown  in  figure  2.  Magnetic  energy  stored  in  the  generator  between  the  coil  and 
armature  is  continually  being  lost  through  the  vacuum/conductor  interface  to  be 
stored  as  magnetic  energy  in  the  conductor  or  dissipated  as  joule  heat.  This  loss 
through  the  interface  can  be  represented  by  effective  resistances  in  the 
equivalent  circuit  (appendix  A). 
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The  circuit  relation  is 
’  (LI)  +  RI  =  0, 


at 

with  solution  at  time  t 


L  and  R  refer  to  total  circuit  values  and  t  =  0  refers  to  the  priming 


condition. 


2.1 


Generator  inductance 


This  is  .a  circuit  representation  for  the  energy  stored  between  the 
armature  and  coil.  Standard  techniques  for  calculating  coil  self  and  mutual 
inductances  are  given  in  reference  (2).  The  geometry  of  figure  1  can  be 
represented  by  a  coil  of  given  length  and  effective  radius,  the  latter  being  based 
on  the  mean  cross-sectional  area  between  the  armature  and  coil.  Alternatively, 
the  coil  can  be  considered  as  a  number  of  interacting  sections.  Further,  the 
corresponding  armature  sections  represent  mirror  images  of  the  coil  sections,  ie, 
the  coil  sections  and  armature  sections  have  similar  ampere  X  turn  values.  In  this 
case  we  have 

N  N  N  r  V 

"  I  I  ^ik  °ik’ 

®  i=i 

where  the  self -inductances  of  the  coil  sections  are 
L.  i  ^  1,  2  . N, 

the  mutual  inductances  between  coil  sections  are 

ifk  =  Ij  2  ......  N, 

and  the  mutual  inductances  between  coil/armature  sections  are 

i  =  1,  2  . N  k  =  1,  2 . M. 

2.2  Armature  effective  resistance 


Because  skin  depths  are  in  general  much  smaller  than  armature  radial 
dimensions  for  the  systems  under  consideration,  a  one-dimensional  field  diffusion 
problem  is  appropriate.  The  expression  for  the  effective  resistance  is  given  in 
appendix  A  as 
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V  -  n  r  „9H 

a  "  -  J  ....(1) 

s 

where  n  is  the  coil  winding  density, 

is  the  magnetic  field  at  the  armature  surface, 
p  is  the  material  resistivity, 

X  is  the  distance  measured  into  the  conductor, 

S  is  the  armature  surface  area. 

Hence,  is  evaluated  from  a  knowledge  of  the  surface  field  and  its 
gradient  and  the  surface  temperature  which  defines  p.  The  surface  field  is  given 
by 

Hs  =  nl. 

During  the  earlier  part  of  the  operation  the  armature  consists  of  two 
sections  having  differing  physical  conditions,  ie,  shocked  and  pre-shocked 
conductors.  In  general,  significant  conductor  heating  occurs  through  joule 
dissipation  and  the  problem  is  non-linear. 

Three  methods  are  employed  to  evaluate  Two  are  based  on 
exponential  approximations  for  the  current  ^and  f'eld)  tin.e  profiles.  The  third 
represents  a  finite  difference  technique  for  arbitrary  surface  profile. 

2.2.1  Equivalent  skin  depth  method 

The  approach  is  similar  to  that  of  reference  (3). 


The  solution  of  the  linear  field  diffusion  equation 


Ug3H/3t  =  p3^H/3x^ 


is  for  the  case  of  an  exponential  applied  surface  field  profile  of  form 

H  =  H  exp(cit) 
s  o 

with  boundary  condition  H(d,t)  =  0,  where  d  is  the  conductor  thickness. 
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H(x,t)  =  (a  exp(-  x/a)  +  (1  -  A)  exp  (x/a)) 
once  the  initial  field  profile  H(x,o)  has  been  overriden. 


....(2) 


Where  a  =  ( - )  is  the  exponential  skin  depth, 


is  the  conductor  permeability. 


A  =  (l  -  exp(-  2d/a)) 


Substitution  of  equation  (2)  in  equation  (1)  gives 
R  =  -nf.p.(2A  -  .1)  s, 

SL  3 


The  joule  deposition  in  the  conductor  is  given  by 


W  =  I 


where  and  k'  =  (  2A  -  1  -  4A(1  -  A)^/a] , 


We  now  make  the  approximation  (3)  that  for  a  d  the  mean  effective 
temperature  rise  is  evaluated  by  supposing  that  2/3  of  the  joule  heat  goes  into  a 
layer  of  thickness  a, 

Sa5E^  =  2/3W.  fit, 
fiEy  =  fiT, 

where  E  is  the  energy  density/unit  volume. 


C  is  the  heat  capacity/unit  volume, 


T  is  the  temperature. 
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For  a  >d  the  deposition  is  taken  as  being  uniform  over  the  conductor  thickness  d, 

Sd6E  =  W.6t. 

V  j 

In  general,  the  exponential  factor  a  is  time  dependent.  Local  values 

defined  by 


a(t) 


1  3h 
H  3t 


are  used  in  the  computations. 


2.2.2  Asymptotic  analytic  method 

This  approach  is  based  on  the  work  of  reference  (4).  For  semi-infinite 
conductors,  analytic  solutions  of  the  non-linear  diffusion  equation 

1J„3H/St  .  |j(p3H/a«) 

can  be  derived  for  the  case  of  very  low  and  very  high  fields  for  particular  time 
profiles  of  the  surface  field.  The  solutions  are  only  valid  near  the  conductor 
surface;  however,  this  is  sufficient  to  satisfy  the  requirements  of  equation  (1). 
Further,  the  solutions  are  only  valid  for  materials  with  properties 

Cy  independent  of  E^, 

....(4) 

p  =  Pg(l  +  kE^). 

These  relations  are  reasonably  consistent  with  real  material  proper¬ 
ties. 


For  a  profile  of  form  exp(ot),  we  have:- 
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(1)  H  «  Hj,, 


H(x,t)  =  H^expC-  x/a), 


u. 


H. 


EyCOjt)  =  -^g. 


....(5) 

....(6) 


(2)  H  » 


H(x,t)=H3(l4|-iJ|^),  ||(x.t)| 

8 


H. 


...(7) 


2 


H  is  the  critical  field  defined  as  (• — y, 


,...(8) 


Ey  is  zero  for  H  =  0, 


Substitution  of  equations  (4),  (5)  and  (6)  or  (4),  (7)  and  (8)  into  (1) 

gives:- 

H  «  Rg  =  Pod  +  H*/Hj)SnVa  =  PofiSnVa. 

H  »  =  Po(H^/H3  +  I  Hg/H^)Sn2/a  =  PofzSnVa. 

The  functions  fi  and  fj  are  shown  in  figure  3.  This  analysis  does  not 
give  values  of  for  intermediate  fields.  However,  figure  3  indicates  that  a 
sensible  extrapolation  may  be  made. 


2.2.3 


Finite  difference  method 


General  solutions  of  the  non-linear  field  diffusion  equation  (3)  and  the 
heat  equation 

3E^/9t  =  pj^  =  p(3h/3x)^, 

3E^/3t  =  C^, 
j  =  current  density, 


can  be  obtained  for  any  surface  profile  Hg(t)  and  any  p(T),  C^vT).  The  appropriate 
finite  difference  equations  are;- 


=  H  +  - ^ - (p  -  p  )(H  -  H  ) 


4p^(6x) 

o 


2  m+i  m-i' 


P  6t 

+  -  (H  . ,  -  2H  +  H  ,  ) , 

(5^)2  n.  m-i 

o 


=  E  + 
vra  vm 


p  6t 
^  *1* 


4(6x) 


(H 

2  in+ 1  m—  1 


The  conductor  is  divided  up  into  a  number  of  equal  increments  5x(n) 

where  the  dependent  variables  have  values  H  ,  E  ,  n  at  time  t 

m'  vm’  ’^m 

X  =  rafix,  d  =  nfix, 

m  =  0,  1,  2 . n. 

The  superscript  +  denotes  values  at  t  +  5t. 


A  reasonable  approximation  to  the  real  situation  is  represented  by  the 
boundary  condition 


H 


n 


=  0. 


t  t 


Coil  effective  resistance 


The  coil  presents  a  far  more  complicated  surface  problem  than  the 
armature.  Let  us  consider  the  two  separate  components  making  up  the  coil 
effective  resistance. 

=  Rj  .  R„, 

where  is  the  component  representing  joule  heating, 

Rj^  is  the  component  representing  magnetic  energy  storage. 

The  calculation  of  R^  for  the  case  of  steady  ac  has  been  the  subject 
of  a  number  of  studies.  The  most  general  case  for  the  linear  problem  has  been 
considered  in  reference  (5)  which  gives  the  relation 

Rj/Rg  =  V  +  n^x(3ui  +  iru2)» 

where  Rg  is  the  ac  resistance  of  the  coil  wire  when  straight,  ie,  when  subjected 
to  skin  effects  alone. 


Vj  Ti>  X>  3j  Y>  '^1  ®^nd  U2  are  functions  of  the  coil  dimensions,  the 
conductivity  and  the  frequency  f  and  account  for  proximity  effects. 

In  general,  one  can  equate  the  resistance  of  a  circuit  carrying  an 
exponential  current  profile  to  an  ac  resistance  of  that  circuit, 

Rj(a)  =  Rj(f  =  a/ir). 

In  order  to  estimate  temperature  rises  and  hence  their  effect  on  R^ 
the  approximation  is  made  that  the  current  flow  is  uniform  in  an  area  of  cross- 
section  such  that 


where  r,„  is  the  wire  radius, 

Rj^  is  the  dc  resistance  of  the  coil. 


Hence,  the  energy  balance  relation  is 

where  %  is  the  coil  wire  length. 

Two  approximations  are  used  to  estimate 

(a)  It  is  well  known  that  for  large  magnetic  fields  and  large  rates  of 
field  change  the  joule  and  magnetic  energy  densities  within  a 
conductor  are  approximately  equal.  Hence, 

(b)  The  coil  is  approximated  by  a  conducting  layer  of  thickness  2r^. 
For  such  a  layer  with  the  boundary  condition  of  H  =  0  at  the  outer 
surface,  the  results  of  section  2.2.1  may  be  applied, 

Sf  -  \  -  k"' ' 

=  ((2A  -  1)  +  4A(1  -  A)2rya), 

S  =  2TTr^i  , 
c  o  ’ 

where  r„  is  the  coil  radius, 
c 

a  is  the  coil  length, 
o 

3.  BASIC  INPUTS  TO  THE  MODEL 


The  electrical  equivalent  circuit  model  of  section  2  will  require  input 
information  over  and  above  the  static  dimensional  description  of  the  generator. 

3.1  Conductor  properties 

The  conductor  property  of  paramount  interest  is  the  electrical 
resistivity.  Because  a  wide  range  of  temperatures  may  be  encountered,  the  data 
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must  cover  many  thermodynamic  states,  ie,  solid,  liquid,  vapour,  plasma.  The 
resistivity  variation  employed  with  the  current  model  is  shown  in  figure  4  where 
the  two  conductors  Cu  and  A1  are  taken  as  examples. 

Reference  (6)  provides  data  in  the  solid  regioii  and  for  a  limited  range 
in  the  molten  region.  The  molten  data  are  linearly  extrapolated  to  the  boiling 
point.  A  further  linear  approximation  is  taken  from  the  boiling  point  to  a 
temperature  of  2  x  10  where  the  plasma  resistivity  model  of  reference  (7)  is 
employed.  Reference  (8)  effectively  gives  the  plasma  conductivity  as 


o  e 

m^e^Z 

e 


SttZ^  ,1, 

- T^(kT)®l% 

(2m  )2Ze^ 
e 


(degenerate  )  (classical  ) 

(plasma  term)  (plasma  term) 

which,  on  inserting  values  for  physical  constants,  reduces  to  (Z  =  1), 

i  =  +  2.92  X  10"®  (Sim), 

P  1.12  X  lO” 

-3\ 

where  is  the  free  electron  density  (m  ). 


Further, 


where  n  .  is  the  density  of  neutrals  +  ions  and  k  is  the  Boltzmann's  constant, 
n+i  ^ 

Figure  4  is  based  on  the  approximation  of  one  free  electron  per  atom. 


Magneto-resistance  effects  are  expected  to  be  small. 


Piezo-resistande  effects  arising  through  magnetic  forces  exerted  on 
the  conductors  are  expected  to  take  the  Bridgnian  form 

_  /V  vy 

P  ~  ’ 

^vo  o 

where  v  and  v^  are  the  specific  volumes  at  the  test  and  zero  pressures 
respectively.  Exponent  values  for  Cu  and  A1  are  2.7  and  3.3  respectively. 
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and  Al. 


Table  1  shows  the  enthalpy  relations  up  to  the  boiling  point  for  Cu 


TABLE  1 


Parameter 

Cu  Value 

Al  Value 

C,3(0) 

3.4  MJ  m~^  deg~* 

2.4  MJ  m  ^  deg  ^ 

®vs 

2.36  X  lo"'*  deg~^ 

5.0  X  lO"**  deg~^ 

C,JMP) 

4.07  MJ  m"’  deg~^ 

2.55  MJ  m  ^  deg  ^ 

®VI. 

0.0 

0.0 

1.6  GJ/m^ 

1.08  GJ/m® 

48  GJ/m’ 

30.5  GJ/m® 

where  the  specific  heat/unit  volume  relations  have  the  form 
^VS  ~  SyjT)  for  the  solid  conductor, 

^VL  "  liquid  conductor, 

T  =  °C, 

^MP  "  point, 

Lp  is  the  latent  heat  of  fusion, 

Ly  is  the  latent  heat  of  vaporisation. 

3.2  Armature  dynamics 

The  time  dependent  geometry  will  be  determined  to  a  large  degree  by 
the  motion  imparted  to  the  armature  by  its  explosive  loading.  The  model 
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described  in  this  report  does  not  generate  this  motion  but  uses  it  as  input.  Such 
information  can  be  generated  either  through  experiment  or  hydrocode  simulation. 
In  the  former  case,  the  expansion  characteristics  are  determined  in  the  main 
from  high  speed  photographic  diagnostics.  Both  framing  camera  photographs  of 
the  armature  and  streak  camera  records  of  suitably  arranged  flashgaps  are  used 
to  build  up  a  description  of  the  funnelling  action.  Such  experiments  are  also 
important  in  determining  the  condition  of  the  conductor  during  its  expansion,  ie, 
how  far  can  it  expand  before  deterioration,  such  as  cracks  which  can  impede 
current  circulation,  sets  in?  Alternatively,  the  expansion  of  the  shock  loaded 
armature  can  be  determined  from  hydrocode  study.  Such  an  approach  will  also 
give  information  on  the  shock  induced  temperatures  in  the  conductors.  Both 
approaches  are  compatible  as  far  as  the  expansion  characteristics  are  concerned. 

The  end  result  of  both  approaches  is  basically  the  cone  expansion 
angle.  In  general,  the  armature  expands  in  the  form  of  a  cone  with  fixed  angle 
whose  truncated  apex  is  displaced  from  the  original  position  of  charge  initiation 
by  the  product  of  detonation  velocity  x  time  from  initiation,  as  indicated  in 
figure  1.  Such  a  characteristic  is  easily  inputted  into  the  model  described  in  this 
report.  In  practice  this  simple  description  is  not  completely  representative.  For 
instance,  expansion  characteristics  will  be  perturbed  by  the  effects  at  the 
armature  ends.  However,  it  is  customary  to  employ  armatures  which  overlap  the 
coil  so  that  such  effects  are  effectively  removed.  Although  the  cone  angle  does 
vary  along  the  cone  at  any  particular  instant  in  time  a  mean  angle  approximation 
represents  a  good  compromise. 

The  explosively  induced  expansion  of  the  armature  represents  the 
only  noticeable  contribution  to  motion  when  magnetic  field  strengths  are 
reasonably  small.  However,  as  fields  increase  electromagnetic  forces  become 
sufficiently  intense  so  as  to  perturb  the  motion  of  the  armature  and  to  induce 
movement  in  the  coil.  Such  forces  play  an  increasingly  important  role  as  the 
generator  is  scaled  up  in  size.  The  model  adopts  the  following  approach. 


Magnetic  fields  are  assumed  to  be  parallel  to  the  conductor  surfaces 
and  produce  forces  normal  to  those  surfaces.  The  coil  is  given  a  field  induced 
expansion  velocity  and  the  armature's  explosively  induced  velocity  has  a  field 
induced  reduction.  Coil  distortion  along  its  axis  is  not  considered.  Because  the 
surface  magnetic  pressure 

Pm  = 

rises  very  rapidly  during  the  period  where  it  is  significant,  then  a  Hugoniot 
approximation  can  be  used  to  relate  particle  velocity  to  the  magnetic  pressure 

Pm  "  Po^s"’ 

where  is  the  density  at  zero  pressure, 

Dg  is  the  shock  velocity, 

u  is  the  particle  velocity. 

Most  metals  have  a  linear  relation  of  the  form 
Dg  =  +  mu. 


where  D  and  m  are  constants, 
o 


The  simple  linear  Hugoniot  relations  for  Cu  and  A1  have  the 
coefficients  shown  in  table  2. 

TABLE  2 


D  » 

o’ 

Po’ 

m/s 

m 

kg/m’ 

Cu 

4000 

1.5 

8900 

A1 

5200 

1.4 

2800 
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4. 


CODE  RESULTS 


The  generator  output  will  depend  on  a  large  number  of  parameters. 
Code  studies  will,  in  general,  be  directed  towards  identifying  the  important 
parameters  and  deriving  optimum  values.  In  this  section  some  typical  examples 
of  output  are  given,  together  with  some  selected  parameter  trends.  A  selected 
comparison  with  experimental  data  is  also  presented. 

A  generator  having  the  following  parameter  description  is  taken  as 
the  basis  for  the  examples:- 

Length:  31.8  mm 

Armature  (Cu)s  radius  7.95  mm;  thickness  0.5  mm;  expansion 

half-angle  14° 

Coil  (Cu):  radius  15.2  mm;  wire  diameter  0.376  mm  (28 

SWG);  pitch  0.442  mm 

Explosive  detonation  velocity:  8800  m/s 

Load;  R  50  mfl,  L  13  nH 

Room  temperature:  20°C;  Shock  release  temperature  of  armature 
170°C 

Priming  current:  200  A 

Unless  otherwise  stated,  the  following  examples  are  directly  relevant 
to  the  above  generator. 

4.1  Typical  results 

Some  of  the  iypical  code  outputs  are  shown  in  figures  5  to  15.  In  this 
example,  Ihe  generator  is  overdriven  with  a  high  priming  current  of  2500  A  in 
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order  to  exaggerate  some  of  the  thermal  effects.  In  each  figure  a  family  of  up 
to  15  curves  are  shown,  each  one  representing  a  different  winding  pitch  value. 
The  time  in  each  figure  is  referenced  to  armature  first  movement  at  the 
beginning  of  the  coil  (it  is  assumed  that  the  generator  is  short-circuited  from  the 
priming  supply  at  this  stage). 

Figure  5  shows  the  number  of  turns  left  in  the  coil.  First  crowbar  of 
the  armature  on  the  coil  is  indicated  by  the  onset  of  the  decrease  in  the  number 
of  turns.  At  6.2  ys  the  generator  has  been  completely  short-circuited. 

Figure  6  shows  the  corresponding  generator  inductance  as  defined  in 

section  2.1. 

Figures  7  and  8  show  the  mean  effective  armature  and  coil 
temperatures  as  defined  by  the  models  in  sections  2.2.1  (equivalent  skin  depth) 
and  2.3.  The  profiles  for  the  armature  in  this  example  have  a  temperature  level 
of  170°C  which  is  derived  from  the  shock  heating  of  the  conductors  from  a  room 
temperature  of  20°C.  Until  ~  3  ys  there  will  be  regions  of  shocked  and  unshocked 
armature;  thereafter  all  the  armature  will  have  been  shocked.  The  profile  in 
figure  7  will  represent  the  hottest  section  of  the  armature.  For  this  example  the 
current  and  hence  the  heating  will  peak  in  the  region  5.5  to  6  ys  so  that  the  steep 
profiles  turn  over  and  show  a  plateau  at  this  point.  Although  the  curves  extend  to 
about  7  ys,  they  have  no  significance  after  6.2  ys  because  the  generator  is 
completely  short-circuited  at  this  time.  Curves  1,  2,  3  and  4  are  different  in 
that  1,  2  and  3  pass  through  the  melting  transition  of  Cu  at  1083°C  and  continue 
into  the  molten  region,  and  4  reaches  the  transition  but  does  not  emerge  from  it. 
The  remaining  curves  are  all  contained  within  the  solid  regime.  Figure  9  shows 
the  latent  heat  curve  that  complements  figure  7.  Curves  1,  2  and  3  acquire  the 
latent  heat  of  fusion  for  Cu  (1.6  G3/m^)  before  heating  terminates,  whereas 
curve  4  only  acquires  0.5  G3/m*  before  heating  ceases.  Figure  8  for  the  coil  is 
similar  to  figure  7  apart  from  the  room  temperature  (20°C)  starting  value,  and 
the  higher  temperatures  reached.  Figure  8  and  its  latent  heat  complement  figure 
10  indicates  that  curves  1  to  12  pass  through  the  melting  transition  and  curve  13 
accumulates  ~  0.5  G3/m*  of  latent  heat  before  heating  stops. 
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For  the  alternative  armature  models  described  in  sections  2.2.2  and 
2.2.3  the  armature  thermal  condition  will  be  described  by  thermal  energy  density 
curves  for  the  surface. 

Figure  1 1  shows  the  armature  effective  resistance  as  defined  by 
section  2.2.1.  It  increases  up  to  first  crowbar  tioiC  through  the  increasing 
shocked  area  of  the  armature  and  the  increasing  exponential  factor  associated 
with  the  current  increase.  After  first  crowbar,  the  resistance  decreases  due  to 
the  reduction  in  armature  area.  However,  because  of  the  extreme  joule  heating, 
curves  1,  2  and  3  show  terminal  resistance  spikes. 

Figure  12  shows  the  joule  component  of  the  coil  resistance  (in  this 
example  the  approximation  R|^;j  has  been  used).  The  profiles  are  similar  to 
armature  resistance,  however,  in  this  case  the  increase  up  until  first  crowbar 
time  is  due  entirely  to  the  increase  in  the  current  exponential  factor.  The  more 
severe  heating  in  the  coil  causes  most  of  the  curves  to  exhibit  the  terminal 
resistance  spike 

The  armature's  magnetically-induced  particle  velocity  is  shown  in 
figure  13.  Because  the  coil  and  armature  conductors  are  both  Cu  in  this  example 
both  will  exhibit  identical  profiles.  The  armature's  magnetically-induced  velocity 
will  be  in  the  opposite  sense  to  the  explosively-induced  particle  velocity  and 
hence  will,  in  addition  to  the  coil  expansion  velocity,  increase  the  time  required 
to  completely  short-circuit  the  generator.  This  point  is  amplified  in  figure  5 
where  it  is  seen  that  the  curves  reach  zero  at  slightly  differing  times.  The 
highest  particle  velocity  in  figure  13  for  curve  1  represents  only  6%  of  the 
explosively-induced  velocity.  In  general,  appreciable  slowing  down  of  device 
operation  only  occurs  for  generators  which  are  significantly  larger  than  the 
example  under  consideration.  Small  generators  are  unable  to  transduce  large 
ff'actions  of  armature  kinetic  energy  into  magnetic  energy  because  of  the  large 
joule  losses  incurred  at  the  high  temperatures  produced. 

Figure  14  shows  the  generator's  current  gain  as  a  function  of  time. 
The  values  are  somewhat  on  the  low  side  because  the  generator  is  overdriven 
with  the  2500  A  priming  current.  The  time  of  peak  varies,  ie,  the  larger  the  gain 
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the  sooner  severe  heating  sets  in  resulting  in  a  turn-over.  Optimum  output  is 
obtained  in  this  example  for  a  wire  diameter:pitch  ratio  of  0:6  (although  the 
optimum  is  not  sharp).  In  general,  this  optimum  occurs  for  close-wound  coils  but 
it  does  tend  to  wander  from  this  position  slightly  when  conductors  get  hot.  It  will 
also  be  noted  that  current  peak  occurs  before  the  generator  is  completely  short- 
circuited.  Although  joule  heating  is  a  significant  factor  in  producing  this  effect 
in  this  example,  it  is  a  general  result  even  for  cold  conductors.  Near  the  end  of 
generator  operation  the  dl/dt  driving  function  (slope  of  figure  6)  is  no  longer 
capable  of  offsetting  the  circuit  resistance  in  the  relation 


dl 


Finally,  figure  15  shows  the  typical  increase  with  time  of  the  current 
exponential  factor  oi(t).  In  this  model  the  factor  retains  its  peak  value  after  peak 
current. 


4.2  Compatibility  of  models 

In  section  2  a  number  of  alternative  models  for  computing  circuit 
parameters  were  given.  In  general,  the  alternative  approaches  are  quite 
compatible.  The  alternative  approaches  for  calculating  generator  inductance  are 
consistent  to  within  ~  1%.  The  first  of  the  two  models  for  calculating  Rj  in 
section  2.3  gives  a  slightly  higher  result.  However,  unless  the  load  impedances 
are  small,  the  current  amplifications  are  not  too  dissimilar. 

An  illustration  of  the  degree  of  agreement  among  the  three  models 
for  computing  armature  effective  resistance  is  shown  in  figure  16  for  the 
example  generator  already  described.  The  peak  current  amplifications  are  shown 
as  a  function  of  priming  current  for  a  wire  diameter/winding  pitch  =  0.8.  Figure 
16  indicates  that  the  equivalent  skin  depth  and  finite  difference  models  are 
virtually  equivalent.  Further,  the  asymptotic  approach  is  consistent  with  the 
other  two  models  at  the  extremes  of  the  range  but  shows  a  divergence  of  up  to 
5%  in  the  intermediate  region.  This  result  is  not  unexpected  and  indeed  the 
divergence  can  be  reduced  by  using  a  more  elaborate  extrapolation  approach  for 
the  asymptotic  model.  The  asymptotic  approach  is  only  correct  in  principle  for 
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H  »  H  »  In  figure  16  the  priming  current  which  leads  to  at  peak 

amplification  is  shown.  The  peak  current  amplification  is  most  strongly 
influenced  by  the  armature  resistance  near  peak  time  where  the  circuit's  time 
constant  is  small.  Hence,  the  asymptotic  approach  is  only  expected  to  be  correct 
for  priming  currents  significantly  above  and  below  ~  200  A  for  this  example.  In 
figure  16  a  linear  extrapolation  for  the  armature  resistance  was  taken  for  the 
intermediate  region. 


f  = 


e 


eH  -  — 
c  e 


{f2(eH^)  -  fi(H^/e)}, 


with  e  =  2.5.  However,  the  result  is  not  very  sensitive  to  the  vaiue  assigned  to  e. 
Although  a  better,  more  elaborate  extrapolation  is  possible,  the  simple  linear 
approach  is  deemed  to  be  satisfactory. 

The  more  exact  finite  difference  approach  is  regarded  as  the 
standard.  However,  because  it  requires  significantly  more  computer  time  it  is 
used  as  a  check,  whereas  the  other  models  are  employed  for  large-scale 
parameter  variation  studies. 


4.3 


Load  effects 


The  example  generator's  performance,  as  a  function  of  load  resis¬ 
tance  and  inductance,  is  shown  in  figure  17.  The  influence  of  load  changes  is  only 
highly  significant  for  low  load  values.  The  predictions  are  compared  with 
experimental  data  for  this  generator  (9).  Agreement  is  seen  to  be  good  over  a 
wide  range  of  load  values  both  in  absolute  value  and  trend. 


4.4 


Influence  of  wire  diameter  and  winding  pitch 


The  generator's  performance  is  maximised  for  an  optimum  combina¬ 
tion  of  wire  diameter  and  winding  pitch.  The  optimising  conditions  for  current 
and  energy  gains  are  in  general  different.  This  effect  is  illustrated  for  our 
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example  generator  in  figure  1?.  Maximum  current  gain  implies  maximum  rate  of 
current  rise.  On  the  other  hand  maximum  energy  gain  ensures  the  highest  output 
current  for  a  given  system  size  (and  priming  energy). 


Figure  18  shows  th^t  for  each  wire  diameter  there  are  optimum  pitch 
conditions  for  current  and  energy  gain.  The  former  lie  in  the  vicinity  of  close- 
wound  coils  while  the  latter  occur  at  larger  pitches.  Further,  the  optimised  wire 
diameter  is  larger  for  energy  gains.  The  optimising  conditions  depend  on  the 
system  parameters.  For  instance,  as  priming  currents,  with  associated  heating, 
are  increased  optimum  conditions  t'nd  towards  larger  wire  diameters  and 
pitches.  This  effect  has  already  been  noted  in  figure  14  and  section  4.1. 

The  example  shown  in  figure  18  indicates  that  the  conditions  for 
optimum  current  and  energy  gains  are  significantly  different.  However,  a 
compromise  choice  can  be  made  which  gives  good  performance  in  both  aspects. 

4.5  Influence  of  expansion  angle  and  expansion  ratio 


The  expansion  ratio  is  defined  as  coil  radius/armature  radius.  In 
figure  19  the  expansion  ratio  is  varied  through  changes  in  armature  diameter 
alone,  ie,  fixed  overall  generator  size.  The  most  striking  result  is  that  output 
improves  significantly  with  expansion  angle.  This  is  attributable  to  the  reduction 
of  generator  run  time.  However,  the  characteristics  show  far  smaller  effects  for 
variation  of  expansion  ratio  at  fixed  expansion  angle.  In  fact,  the  energy  gains 
show  something  of  an  independc»nce  of  expansion  ratio.  This  result  is  not 
unexpected  because  the  most  significant  transduction  of  energy  and  loss  through 
diffusion  occurs  during  the  latter  stages  of  the  generator's  operation.  During 
these  final  stages,  the  generator's  geometry  will  be  independent  of  the  initial 
expansion  ratio.  Because  in  practice  one  has  to  maintain  the  physical  integrity  of 
the  expanding  armature  in  order  to  maintain  its  good  current  carrying 
characteristics,  a  greater  expansion  ratio  will  be  associated  with  a  thicker 
heavier  conductor  and  less  explosive.  This  would  suggest  a  reduction  of  expansion 
angle  and  hence  some  drop  in  energy  gain  with  increasing  expansion  ratio.  A 
lower  limit  for  the  desirable  expansion  ratio  would  be  that  set  by  the  t.  ondition 
for  the  armature  to  reach  its  terminal  expansion  velocity.  This  limit  is  of  the 


order  of  1.5. 


A  great  deal  of  effort  has  been  expended  in  the  past  in  order  to 
increase  the  expansion  ratio  significantly.  The  above  result  suggests  that 
alternatively  one  should  look  for  minimum  conductor  thickness  at  more  modest 
expansion  ratios. 

4.6  Influence  of  priming  current 

The  CMF  generator  described  in  this  report  is  essentially  a  current 
amplifier.  At  low  priming  current  levels,  the  current  and  energy  amplifications 
are  independent  of  priming  current.  However,  at  higher  levels,  joule  heating 
becomes  significant  leading  to  a  loss  of  performance.  In  fact  the  output  shows  a 
limit  in  total  current.  However,  by  this  stage  the  generator/priming  assembly  is 
working  inefficiently  on  an  output/volume  basis  and  a  bigger  generator  or  one 
with  a  different  coil  is  required.  At  the  higher  current  levels,  magnetic  pressure 
may  also  contribute  to  the  loss  in  amplification.  This  influence  will  be  small  in 
small  generators.  Some  of  these  features  are  illustrated  in  figure  20.  When  the 
total  output  current  for  the  28  SWG  coil  has  reached  its  upper  limit,  the 
generator's  current  amplification  has  dropped  to  about  a  third  of  its  low  current 
value.  Figure  20  indicates  that  for  this  small  generator  the  magnetic  pressure 
contribution  to  the  loss  of  amplification  is  small.  Additionally,  figure  20 
illustrates  how  optimum  wire  diameter  increases  with  priming  current. 

Good  agreement  is  shown  between  the  predictions  and  experimental 

data  (9). 

4.7  Influence  of  resistivity 

One  of  tne  material  physical  properties  that  has  a  significant 
influence  on  generator  output  is  the  conductor  resistivity.  This  applies  both  to 
the  coil  and  armature. 

Figure  21  shows  the  generator  output  as  a  function  of  priming  current 
for  various  ambient  temperatures,  both  above  and  below  normal  room  tempera¬ 
ture.  The  curves,  which  show  a  significant  temperature  dependence  at  a  low 
priming  level,  converge  for  higher  priming  currents  where  the  initial  tempera¬ 
tures  differences  are  offset  by  joule  heating.  However,  in  practice,  changes  in 
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ambient  temperature  are  likely  to  produce  effects  in  addition  to  resistivity 
changes.  For  instance,  explosive  performance  may  change  and  differential 
expansions  will  produce  a  poorer  fit  between  the  armature  and  its  explosive 
charge,  especially  at  very  low  temperatures. 

Figure  22  shows  resistivity  effects  through  the  use  of  different 
conductor  materials  in  the  coil  and  armature.  All  four  combinations  for  the  two 
conductors  Cu  and  A1  (armature  and  coil)  are  taken  as  examples.  The  results,  for 
two  values  of  priming  current,  show  the  expected  significant  loss  of  output  as 
one  moves  from  Cu  to  Al.  However,  the  results  for  the  A1  armatures  will  be 
a  little  improved  in  practice  because,  being  somewhat  lighter,  they  will  have 
slightly  larger  expansion  angles.  The  results  in  figure  22  are  all  for  an  expansion 
angle  of  14°. 

5.  CONCLUSIONS 

A  model  for  a  compressed  magnetic  field  generator  employing  a 
round  wire  single  pitch  helical  coil  has  been  presented.  Magnetic  diffusion  losses, 
including  joule  heating,  in  both  armature  and  coil  are  taken  into  account, 
together  with  the  dynamic  effects  of  magnetic  pressure. 

The  aim  of  the  computational  side  has  been  to  avoid  undue 
mathematical  complexity  while  still  incorporating  the  essential  physical 
features.  Although  alternative  techniques  are  described  for  computing  circuit 
parameters,  these  lead  to  consistent  results.  As  a  general  guide,  ten  complete 
runs  on  an  IBM  370/168  machine  take  ~  1  min.  However,  the  duration  may  be 
increased  by  an  order  of  magnitude  when  the  finite  difference  model  is  employed 
for  computing  armature  effective  resistance.  This  more  exact  method  is  usually 
used  only  as  a  "spot  check".  The  modest  run  time  enables  one  to  carry  out 
extensive  parametric  studies  quickly  and  cheaply  with  a  view  to  minimising  the 
number  of  time  consuming  and  costly  experiments. 

The  good  agreement  between  prediction  and  experimental  data  is 
particularly  encouraging. 
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APPENDIX  A 


Consider  a  generator  having  current  I  in  a  coil  of  n  turns/unit  length. 


The  magnetic  energy  stored  in  vacuo  is 

"m  “  (  T 


(Al) 


where  the  field  is  a  function  of  time  and  position  and  the  integral  is  take  over 
the  volume  bounded  by  the  generator  conductors. 


From  the  definition  of  inductance  L, 


Differentiation  of  equations  (Al)  and  (A2)  yields 
®  (H.H)<1V  +  Jg~(H.H)u.d?, 

=  Lli  +  Ili^ 
m  2  ’ 

where  u  is  the  conductor  velocity  and  S  the  bounding  surface. 


Equating  terms  due  to  field  and  volume  changes 

Lli  =  ;  Hi 

V 

=  J  ii(H.H)u.dl. 

^  S 


The  general  circuit  relation  is 
^  (LI)  +  RI  =  0 

and  the  energy  balance  relation 

-  ri2  =»  W  +  i  LiS 

m  4 

=  W  +  J  (H.H)u.dS. 
in  _  2 


(A2) 


(A3) 


(A4) 


Further,  for  the  vacuum  where  electrostatic  energy  storage  is  negligible  and 
where  there  is  no  dissipation  we  have 

J(E  X  H).dS  =  -  A  L(H.H)dV.  (A5) 

S  V 

Substitution  of  equations  (A3)  and  (A5)  into  (A4)  gives 

-  RI*  =  -  J(E  X  H).ls  +  J  ]i  (LH)u.dS.  (A6) 

S  S  ° 

The  appropriate  field  equations  for  the  moving  conductor  are  (using 
standard  nomenclature) 

V  X  H  =  j, 

®EFF  =  ^  ^  ^  ^  V- 
These  may  be  manipulated  to  give 

E  X  H  =  P(V  X  H)  X  H  +  U^(^.H)u  -  y^(H.u)H, 
which  on  substitution  in  equation  (A6)  yields 

-  RI^  =  -  J  p(V  X  H)  X  H.ts  +  J  p  (H.u)H.’3s,  (A7) 

S  S 

Hence,  the  energy  lost  through  the  vacuum/conductor  interface  to  be 
dissipated  as  joule  heat  or  stored  as  magnetic  energy  in  the  conductor  can  be 
evaluated  through  an  equivalent  circuit  or  effective  resistance  R  of  the  form 
depicted  in  equation  (A7). 

If  H  and  u  or  H  and  dS  are  everywhere  perpendicular,  then  the  second 
term  in  equation  (A7)  is  zero. 

Further,  for  an  infinitely  long  system 
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JOULE  COMPONENT  OF  COIL  RESISTANCE  Rj 
AS  FUNCTION  OF  TIME 
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OUTPUT  AS  A  FUNCTION  OF 
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DIFFERING  AMBIENT  TEMPERATURES 
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DIFFERING  CONDUCTOR  MATERIALS 
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Quantity 

Unit 

Symbol 

Conversion 

Basic  Units 

Lenxjih 

metre 

■ 

1  m  -  3,2808  ft 

1  ft  -  0.3048  m 

Mass 

kilogram 

kg 

1  kg  -  2.204A  lb 

1  lb  -  0.45359237  kg 

1  ton  “  1016.05  kg 

Derived  Units 

Force 

newton 

N  -  kg  m/s^ 

1  N  -  0,2248  Ibf 

1  Ibf  -  4.44822  N 

Work,  Enerjty,  Quantity  of  Heat 

joule 

J  -  N  m 

1  J  -  0.737562  ft  Ibf 

1  J  -  9,47817  »  lO-'*  Btu 

1  J  -  2.38846  »  lO""*  kcal 

1  ft  Ibf  -  1,35582  J 

1  Btu  -  1055.06  J 

1  kcal  -  4186.8  J 

Power 

watt 

W  -  J/s 

1  W  -  0.238846  cal/a 

1  cal/a  -  4.1868  U 

Electric  ChatRe 

coulomb 

C  •  A  s 

- 

Electric  Potential 

volt 

V  -  W/A  -  J/C 

- 

Electrical  Capacitance 

farad 

F  -  A  s/V  -  C/V 

- 

Electric  Resistance 

ohm 

n  -  V/A 
s  -  1 

- 

(inductance 

seamen 

- 

Magnetic  Flux 

waber 

Wb  •  V  a 

- 

Magnetic  Flux  Density 

tesla 

T  -  Wb/m* 

- 

Inductance 

henry 

H  -  V  s/A  -  Wh/A 

• 

Conolex  Derived  Units 

Angular  Velocity 

radian  per  second 

rad/s 

1  rad/e  -  0.159155  rav/s 

1  rav/s  “  6.28319  rsd/s 

Acceleration 

metre  per  square  second 

m/s* 

1  m/s*  -  3.28084  ft/s* 

1  ft/s*  •  0,3048  m/o* 

Angular  Acceleration 

radian  per  square  second 

rad/s* 

- 

Pressure 

neuton  per  square  metre 

N/m^  -  Pa 

1  N/m*  -  145.038  -  10"^  Ibf/ln* 

1  Ibf/in?  -  6.89476  »  10*  N/m* 

bar 

bar  -  10*  N/m^ 

1  In.  Hr  -  3386.39  N/m* 

Torque 

newton  metre 

V  m 

1  N  m  -  0.737562  Ibf  ft 

1  Ibf  ft  -  1.35582  N  m 

Surface  Tension 

newtor.  per  metre 

N/m 

1  N/m  -  0.0685  Ibf/ft 

1  Ibf /ft  -  14.5939  N/m 

Dynamic  Viscosity 

newton  second  per  square  metre 

N  s/m* 

1  N  s/m*  -  0.0208854  Ihf  s/ft* 

1  Ibl  s/ft*  -  47.8803  S  a/m* 

Kinematic  Viscosity 

square  metre  per  second 

m*/8 

1  m*/B  -  10.7630  ft*/B 

1  ft*/s  -  0,0929  m*/B 

Thermal  Conductivity 

watt  per  metre  kelvln 

W/m  K 

Odd  Units* 

Radioactivity 

hecquarel 

Bq 

1  Bq  -  2.7027  »  10"'>  Cl 

1  Cl  -  3.700  •  10*®  Be 

Absorbed  Dose 

gray 

Cy 

1  r.y  »  100  rad 

1  rad  -  0.01  Cy 

Dose  Equivalent 

slevert 

Sv 

1  Sv  «  lOO  ram 

1  ram  •  0.01  Sv 

Exposure 

coulomb  per  kilogram 

C/kg 

1  C/kg  •  3876  R 

1  R  -  2.58  »  10—  C/kg 

Rata  of  Leak  (Vacuum  Syatama) 

millibar  litre  par  second 

mb  1/s 

1  mb  •  0.750062  torr 

1  torr  -  1.33322  mb 

•t>)eaa  tanaa  ara  raeoanlaad  tarao  within  the  aatrle  ayatm. 


